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We present a theoretical study of diluted magnetic semiconductors that includes spin-orbit cou- 
pling within a realistic host band structure and treats explicitly the effects of disorder due to 
randomly substituted Mn ions. While spin-orbit coupling reduces the spin polarization by mixing 
different spin states in the valence bands, we find that disorder from Mn ions enhances the spin 
polarization due to formation of ferromagnetic impurity clusters and impurity bound states. The 
disorder leads to large effects on the hole carriers which form impurity bands as well as hybridizing 
with the valence band. For Mn doping 0.01 x < 0.04, the system is metallic with a large effective 
mass and low mobility. 

PACS numbers: 71.23.-k, 72.25.-b, 75.10.Lp, 75.50.Pp, 75.30.Hx, 78.30.Ly 



Ferromagnetic semiconductors with high transition 
temperatures and large polarization are important for 
spintronics applications. The III-V diluted magnetic 
semiconductors (DMS) showing fcrromagnetism with 
Curie temperatures Tc ~ 170K^-^ hold considerable 
promised It is generally agreed that fcrromagnetism be- 
tween the magnetic impurities is mediated via the in- 
teraction with itinerant carriers^ however, a comprehen- 
sive theoretical understanding has been hampered by the 
complexity arising from the interplay of electron-local 
moment interaction, disorder from the doped Mn ions, 
and spin-orbit (SO) effects. 

Theoretical predictions based on the density func- 
tional theory (DFT)^i^ii^ have predicted half-metallicity 
of Gai_a;Mna;As and other DMS. These calculations are 
performed with a supercell containing a small number of 
atoms, often with only one Mn ion in the simulation cell, 
ruling out any disorder effects. Other methods, such as 
the KKR-CPA-LDA method^ (Korringa-Kohn-Rostoker 
method with coherent potential approximation and lo- 
cal density approximation of density functional method), 
include disorder effects only within an effective medium 
approximation. In addition, most of DFT calculations 
are performed without including SO coupling and, as 
a result, overestimate the spin-polarization of carriers. 
Other theoretical approaches based on single-band tight- 
binding metho d^^d^d^ can treat disorder effects prop- 
erly but do not include SO couphng effects. Just as 
the DFT method, these approaches tend to overesti- 
mate the spin-polarization. Multiband methods based 
on k • p approximation^'^' includes the SO coupling 
effects and realistic band structure for the host but the 
disorder effects are completely ignored. 



In this Letter, we address the interplay of disorder 
and SO effects on the electronic structure and the spin- 
polarization of DMS by using a realistic multiband tight- 
binding model. The magnetic coupling between the lo- 
calized Mn substitutional impurity electrons and the itin- 
erant host electrons is included via an exchange param- 
eter. We find the hole states have dual characteristics 
with weight in an impurity band as well as in the valence 
band of the host GaAs, resolving the controversial issue of 
the location of the holes donated by Mn. The density of 
states for low impurity concentration (~ 1% Mn) shows 
deep in-gap localized states at Mn sites. Even at low 
doping, strong hybridization with the valence band states 
pulls considerable weight in the gap region which remains 
separated from the impurity band. For higher dopings of 
about 4% there is some mixing between impurity and 
valence band states, but the impurity states still retain 
their distinct characteristics and are more localized than 
the valence band states.. These states around the impu- 
rity band have a larger effective mass and lower mobil- 
ity than parent semiconductor valence band states, and 
have recently been observed through the optical transi- 
tions between valence and impurity states in midinfrared 
regioni^. 

Another very surprising feature of our results is that 
spin polarization is larger in disordered systems com- 
pared with ordered systems for the same Mn concentra- 
tion, due to the formation of strongly polarized bound 
states and Mn clusters. Similar enhancement of the 
Curie temperature was found within a single-band tight- 
binding model for DMSji^ 

Model: We start with the sp-d exchange model H = 
Hq + Hex , containing the Hamiltonian of the host semi- 
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conductor in the clean limit Hq and the effective ex- 
change coupling between the localized and itinerant spins 
iJea;.— The electronic structure of GaAs is described by 
the Effective Bond Orbital Model (EBOM) developed by 
Chang. The effect of the magnetic impurities is de- 
scribed by 



Hp, 



/,n,m 



Cim) 



(1) 



where cr is the Pauli spin matrix, and S/ are the lo- 
cal spin moments arising from the impurities /, which 
we treat as classical spins of magnitude 5=5/2. n and 
m are the combined indices of the spin and the bond 
orbital, i.e., n = {r, a}, where in the current 8-band 
model T^lt,!} and a—{s,Px,PyjPz}- The impurities 
are placed at the EBOM fee lattice sites and the ma- 
trix elements are calculated in the effective bond basis. 
The exchange couplings capture the hybridization of im- 
purity states with itinerant electron states as well as the 
interaction between electrons mediated by impurity sites. 
Given that Mn atoms in GaAs are substitutional impuri- 
ties with site symmetry of the Ga sites, it can be readily 
shown that the exchange coupling does not mix different 

By comparing 

the single impurity bound state energy with band edge 
emission measurementa^^ of 113 meV, we determine the 
exchange coupling between the valence band holes and 
Mn 3d electrons, Jpd = Jtp,tp = -2.48 eV. The ferromag- 
netic exchange coupling of conduction electrons with Mn 



3d electrons, Jsd = Jt 



are known to be very weak. 



generally of order 0.1 eVi^ We find our results are in- 
sensitive to this coupling and hereafter we set Jsd = 0.1 
eV. 

The advantage of our EBOM based sp-d kinetic ex- 
change model is that both the realistic band structure of 
host materials, including the SO effects, and the disorder 
effects are treated accurately and on an equal footing. 
EBOM not only reproduces k ■ p model results at the 
zone center but also yields much better band structure 
away from the zone center. The disorder effects due to 
the random position of Mn ions are captured exactly for 



a given realization of the Mn ion configuration. 

Bound states of Mn pairs: Before we look into many- 
Mn ion configurations, we investigate the carrier states 
of a Mn ion pair, simulated using a supercell contain- 
ing 586 GaAs primitive cells and two Mn ions, shown 
in Fig. [TJ Note that the carrier states of a Mn ion pair 
are very strongly bound, located well inside the GaAs 
band gap region, i.e. energy region above eV in Fig. [T] 
The nearest neighbor Mn impurity pair forms the deepest 
bound state with an energy « 0.55 eV above the GaAs 
valence band maximum. The spin-resolved local den- 
sity of states (DOS) of Mn pairs shows that the bound 
states are completeley spin-polarized. Bound states of 
large MnAs clusters of more than two Mn ions generally 
induce even deeper bound states (not shown here). The 
non-monotonic distance dependence of the pair bound 
states indicates a strong anisotropy of the ferromagnetic 
coupling due to the zinc-blende crystal symmetry. 

We next analyze the effect of a finite concentration of 
impurities. 

Disorder effects: Randomly substituted Mn ions are 
included explicitly by diagonalizing the Hamiltonian H 
and evaluating the doped hole states for a given local 
spin configuration. Properties, such as the local DOS 
and conductivity are calculated for a given Mn configu- 
ration are then averaged over a large number of Mn ion 
configurations. We assume a coUinear magnetic ground 
state where all localized magnetic ion spins are aligned 
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FIG. 1: Binding energy of carriers to Mn-Mn pairs. Mn^ 
indicates the ith nearest-neighbor Mn pair configuration. The 
horizontal dashed line indicates the position of the valence 
band maximum of host GaAs. 



FIG. 2: Density of States (DOS) of Ga-i-^Mn^^ks. (a) and 
(c): EBOM with randomly placed Mn ions, (b) and (d): 
EBOM with one Mn ion in a supercell. (a) and (b) are for 
X = 0.010 and (c) and (d) are for x = 0.037. Solid (red) and 
dashed (blue) curves correspond to the majority and minority 
spin DOS. 100 disorder configurations were used for the dis- 
order average. Energies are relative to the host GaAs valence 
band maximum energy. The dashed vertical lines indicate the 
position of the chemical potential for uncompensated charge 
density.— Notice that magnetic impurities generate a large 
DOS in the GaAs band gap region [0, 1.52] eV due to strong 
hybridization with the valence band. At both dopings there 
is a peak in the DOS at ~ 0.55 eV (indicated by arrows in (a) 
and (c)) originating from the bound states of nearest neighbor 
Mn ion pairs. 
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along the magnetic easy axis along (100) direction. 

Fig. [2] (a) and (c) show that random locations of the 
Mn impurities generate long tails and a peak around the 
bound state energy of Mn pairs (w 0.55 eV) in the GaAs 
gap region [0, 1.52] eV. In supercell calculations with the 
same impurity concentration, Fig.[2](b) and (d), the DOS 
has large weight in the GaAs gap region but does not 
show the long tails or the extra peak, seen in disordered 
cases. The impurity band around 0.55 eV is separated 
from the main bands at low doping concentrations of x = 
0.01. At a; = 0.037, the impurity induced states still show 
a strong peak but get more strongly hybridized with the 
valence bands. 

The deep in-gap states in disordered impurity configu- 
ration originate from the bonding states of neighbor im- 
purity pairs. The position of the largest impurity band 
DOS, indicated by arrows in Fig. [5] (a) and (c), corre- 
sponds to the deepest bound states of nearest neighbor- 
ing Mn pairs in Fig. [U « 0.55 eV. The deeper states at 
^ 0.7 eV in Fig. [5] (c) at a;=0.037 comes from the bound 
states of cluster of more than two Mn ions. In marked 
contrast, calculations based on effective medium theories 
fail to capture the presence of clusters. In coherent poten- 
tial approximation, for example, the impurity-originated 
states give rise to the impurity band spectral weight near 
the single impurity bound states, which is « 110 meV in 
Gai-xMn^^As , but the long tails deep inside the band 
gap region are not found. 



The presence of in-gap states is of great importance 
for transport and optical properties as they influence the 
inter- and intra-band transition. To study the nature of 
the impurity bands, we calculate the inverse participation 
ratio, IPR{e) — J (ir|V'c(r)|^, where e is the eigenen- 
ergy of the Hamiltonian corresponding to the wavefunc- 
tion tpe{^)- For extended states, IPR{e) scales with the 
system size as ~ ^/V, where V is the volume of the sys- 
tem whereas for localized states, IPR{e) is independent 
of the system size. 

Fig. [3] shows that at impurity doping concentrations 
x=0.037 the carrier states at the Fermi energy are ex- 
tended; hole states with energy < 0.25 eV scale nearly 
linearly with the size of the system, indicating that car- 
riers are delocalized. When the carriers are compensated 
as much as 50% by other defects, such as As antisites, the 
wavefunctions are not completely extended. This local- 
ization tendency might explain the heavier hole effective 
mass in recent midinfrared measurements.^^ For a;=0.01, 
on the other hand, the states at Fermi level are not com- 
pletely extended even without charge compensation. For 
both impurity concentration, the deep in-gap states are 
strongly localized. Our IPR study indicates that the 
charge compensation due to other defects, such as Asca 
antisites, could increase the hole effective mass and lower 
the hole mobility if the chemical potential shifted signif- 
icantly. 

Spin-polarization: Considerable experimental effort 
has been invested in achieving higher ferromagnetic tran- 
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FIG. 3: Inverse participation ratio (IPR) of Gai_j,Mna; As . 
Left panels (a) and (c), are for x = 0.010 and right panels 
(b) and (d), are for x = 0.037. (a) and (b): IPR{e) plotted 
in logarithmic scale as a function of the carrier energy, e. (c) 
and (d): scaling of IPRn, where A'^ is the size of the system. 
The horizontal lines show the IPRn2 /IPRni oc N1/N2, ex- 
pected for extended wavefunctions. The energies are relative 
to the GaAs valence band maximum located at eV. The 
dashed vertical lines denotes the position of the chemical po- 
tential at given charge fraction p = (hole density/Mn density). 
Strong size dependence indicates extended states, whereas lo- 
calized states show little size dependence. For each system 
size, the participation ratio was averaged over 100 disorder 
realizations. 
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FIG. 4: Disorder effects on spin polarization, (a) Density 
of states (DOS) of GaAs with the magnetic impurity con- 
centration X = 0.037 and the exchange coupling Jpd — -1.5 
eV. Upper panel: disordered. Lower panel: ordered array. 
Red(solid) and blue (dashed) lines denote majority and mi- 
nority spin DOS, respectively. The energy is relative to the 
undoped GaAs valence band maximum. The vertical lines 
are the position of the chemical potential. The hole states 
with energy above the chemical potential are occupied, (b) 
The carrier spin-polarization as a function of impurity con- 
centration for two values of Jpd = -2.48 eV (triangles) and 
Jpd = -1-5 eV (circles). The filled (empty) symbols denote 
disordered (ordered) systems. 
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sition temperatures and higher carrier spin-polarization. 
We next discuss the dependence of the carrier spin po- 
larization on the material parameters and impurity con- 
centration. 

Our calculation shows that holes in Gai„a;Mna;As are 
nearly half-metallic for Mn impurity concentrations 
0.01 ^ a; < 0.05 due to the strong exchange coupling, 
Jpd= -2.48 eV. In Fig. [2] (a) and (c), deep bound states 
of Mn clusters pull the majority spin DOS far into the 
band gap region which fully polarizes the carriers. The 
spin-resolved DOS of the ordered system, Fig. [2] (b) and 
(d), also shows nearly complete spin-polarization. Thus 
SO coupling has little influence on the carrier spin polar- 
ization for systems with strong exchange coupling. 

The interplay of SO coupling and disorder is more pro- 
nounced in systems where the coupling between the lo- 
calized spin moment and the carrier spins | Jpd \ is weaker 
relative to the SO coupling A. To simulate systems with 
a smaller | Jp^^j/A, e.g., In2,Mni_j.Sb, we use GaAs mate- 
rial parameters, i.e., the same SO coupling, but a smaller 
exchange coupling. Fig. [1] (a) shows the DOS of GaAs 
with the magnetic impurity concentration of a;=0.037 and 
the exchange couphng Jpd = -1.50 eV. Within our model, 
Jpd = -1.50 eV is too weak for a single impurity to in- 
duce bound states but is capable of forming bound states 
of Mn clusters with non-zero spin-polarization. For a 
weaker exchange coupling, SO effects generate consider- 
able weight in the the minority spin DOS, giving a sub- 
stantial contribution to the polarization for all magnetic 
impurity concentrations. We further see, in Fig. 2] (b). 



the relative enhancement of spin-polarization in disor- 
dered systems, compared with the ordered case which 
grows with increasing impurity concentration. For a 
strong exchange coupling, on the other hand, the spin 
polarization reaches almost 100% regardless the impu- 
rity concentration. 

Our results reveal an important aspect of DMS in spin- 
tronics applications. The spin-polarization and mobility 
of the carriers in a DMS are determined not only by the 
exchange-coupling but also by the geometrical arrange- 
ment of magnetic impurities. While a stronger exchange- 
coupling increases the spin-polarization, it also results in 
more tightly bound states of Mn clusters and, eventually, 
leads to low mobility of holes. Our result also indicate 
that excessive amounts of impurities are unnecessary for 
optimal spin-transport applications. With increased im- 
purity concentration, there are more chances of large im- 
purity clusters, which contribute to spin-polarization but 
degrade the carrier mobility by forming deeper bound 
states. To balance the large spin-polarization and the 
high mobility of itinerant carriers, a delicate modula- 
tion of exchange-coupling and impurity randomness is 
required. 
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